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, We show that one loop GUT scale threshold corrections to gauge couplings 
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JL are a significant constraint on the GUT symmetry breaking sector of the the- 
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ory. The one loop threshold corrections relate the prediction for a s (Mz) to 



X 
S3 



the proton lifetime. We have calculated these corrections in a new complete 
SO(10) SUSY GUT. The results are consistent with the low energy measure- 
ment of a s (Mz)- We have also calculated the proton lifetime and branching 
ratios in this model. We show that proton decay rates provide a powerful test 

for theories of fermion masses. 
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I. INTRODUCTION 



It has been known for some time now that Supersymmetric [SUSY] Grand Unified The- 
ories [GUTs] § 

• suppress the non-SUSY GUT prediction for proton decay, via heavy gauge boson 
exchange, by increasing the GUT scale, and 

• lead to a 10% increase (for one pair of Higgs doublets) in the GUT prediction for the 
weak mixing angle, sin 2 6w, compared to non-SUSY GUTs. 

Now, since sin 2 9\y is known to better than 1% accuracy, one naturally uses the values 
of the two low energy parameters, the fine structure constant, a, and sin 2 9w as input to 
determine the GUT coupling, ocg, and the GUT scale, Mq, and predicts the value of the 
strong coupling, a s . [Note, all low energy parameters are typically evaluated at the scale, 
Mz] This prediction is in excellent agreement with the low energy data The global 

fit to a s from measurements at all energies gives a s = 0.117 ± 0.005 with low (high) 
energy measurements favoring values of a s < (>) 0.12. A global fit to all electroweak 
data § gives a s = 0.127 ± 0.005 ± 0.002 ± 0.001. In comparison, SUSY GUTs predict |J 
a s = 0.127±0.005±0.002, where the first error is the uncertainty in the low energy sparticle 
spectrum, and the second, the top and Higgs masses. Note that the non-supersymmetric 
GUT prediction |a s = 0.073 ± 0.001 ± 0.001 is clearly ruled out. 

In addition, SUSY GUTs contain new sources of proton decay (dimension 5 operators) 
mediated by color triplet Higgs fermions |||7|. They typically lead to the dominant decay 
mode p — > K + v with the lifetime, r p oc M t 2 , where M t is an effective color triplet mass. 

These results are valid for any grand unified symmetry group, as long as SU(3)x SU(2)x 
U(l) is embedded in an unbroken simple subgroup above the GUT scale. Consideration of 
fermion masses, however, leads us to the special grand unified symmetry group SO(10) [§. 
It is the unique group which combines all the fermions of one family {u, d, e, is} into one 
irreducible 16 dimensional representation with the addition of only one electroweak singlet 
state, v, the so-called right-handed neutrino. Consequently, SO(10) Clebschs can be used 
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to relate all charged fermion mass matrices. In this context, a paper by Anderson et al. 
0] [Paper I] showed how all charged fermion masses and mixing angles can be reasonably 
described in terms of four dominant effective mass operators at Mq. Hall and one of us 
(S.R.) [Paper II] then showed how to extend the effective theory at Mq to a complete 
SO (10) SUSY GUT valid up to the Planck or string scales, M. That theory included several 
sectors necessary for - (1) GUT symmetry breaking, (2) doublet-triplet splitting (Higgs 
sector), (3) charged fermion masses (fermion mass sector) and also (4) giving a large mass 
to the electroweak singlet neutrinos (neutrino mass sector). It contained sufficient symmetry 
to be "natural." Hence, it included all operators consistent with the symmetries and only 
contained the SO(10) states[] — 

(n i6 + 3) 16 + me 16 + Tiio 10 + n 45 45 + n 54 54 + m 1 (1) 

with n 16 , n 10 , n 45 , ra 54 , specific integers. Finally, the effective mass operators of paper I 
were recovered upon integrating out states with mass greater than Ma- 
in any complete GUT, the prediction for a s receives corrections at one loop from thresh- 
olds at the weak scale, including the sparticle spectrum, and at the GUT scale from states 
with mass of order Mq. The dominant contribution at Mq comes from states in the sym- 
metry breaking and Higgs sectors. 

The Higgs sector, on the other hand, also affects the proton decay rate. The effective 
color triplet Higgs fermion mass, M t , must be significantly larger than Mq in order to 
increase the proton lifetime. The Higgs doublets on the other hand must remain massless 
at the GUT scale. 
In this paper, 

• we show that the threshold corrections due to doublet-triplet splitting, in the simplest 
models, increases the predicted value for a s (M z ). Thus, while the GUT predictions for a s 
and r p are unrelated at the tree level, at one loop the two effects are coupled. One cannot 



1 These two features are consistent with a stringy origin to such a model [11| 



suppress proton decay without at the same moment increasing the prediction for a s . We 
then show that the model of paper II is ruled out by this effect. 

• We also present a new complete SUSY GUT which appears to be consistent with all 
low energy data at the la level ||12|| . This model has in fact a much simpler GUT symmetry 
breaking sector than that in paper II with fewer states. We present some typical results for 
proton decay rates in this model. In a future paper [DJ, we will present a more detailed 
study. 

II. ONE LOOP THRESHOLD CORRECTIONS AT M G . 



One loop threshold corrections at M G for gauge couplings are given by ]TJ 



a- l (M G )=a G l -A l (2) 
where Aj is the leading log threshold correction to a«. Furthermore, 

where the sum is over all superheavy particles and b\ is the contribution the superheavy 
particle ( would make to the beta function coefficient hi if the particle were not integrated 
out at Mq. 

At one loop the definition of the GUT scale is somewhat arbitrary. In order to avoid 
large logarithms it must certainly be at the geometric mean of the heavy masses; otherwise 
we are free to choose its value. A particularly convenient choice is to define M G as the scale 
where the two gauge couplings, a,, i — 1, 2, meet. At this point Ai = A 2 and we define 

a G = ai (M G ) =a 2 (M G ). (4) 

Then define 

e 3 = (a 3 (M G ) - a G )/a G , (5) 

i.e. the relative shift in a 3 at M G . In general, a value of e 3 ~ — (2 — 3%) is needed to obtain 
a, ~ 0.12. 



A. Formula for £3 in a general SO(10) theory 



The most general SO(10) SUSY GUT we consider includes the complete multiplets listed 
in equation (|l|) with arbitrary values of uiq, nio, n^, n 54 , n\. These states can be decomposed 
into their SU(3) x SU(2) x U(l) content by first considering the decomposition to SU(5) 
[see table [| and then the rest of the way [[n| . 

We have derived a general formula for e 3 in this case (see appendix 1) given by the 
expression (valid to lowest order in die) 

J M 1 
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where M' d is defined as the doublet mass matrix with the massless Higgs doublets pro- 
jected out.0 Note, €3 only depends on the number of states in the theory through the mass 
matrices M 7 in each charge sector 7. Moreover, the effective determinants, defined above, 
explicitly take into account the special contributions of vector multiplets to the threshold 
corrections. 

Putting in the values of b] for i = 1, 2, 3 into eqn. (0) we find 
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The dominant contribution to €3 comes from the GUT symmetry breaking and elec- 
troweak Higgs sectors with superspace potential W sym breakin g and Wmggs, respectively. An 



2 Our notation for the states discussed above is found in table ||. 
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additional small contribution comes from the fermion mass sector. These sectors are by 
reasons of "naturalness" necessarily invariant under several U(l) symmetries (including an 
R symmetry). In appendix 2 we show that these symmetries impose stringent constraints 
on the form of €3. In brief, €3 is implicitly a function of the vacuum expectation values 
[vevs] of fields in the symmetry breaking sector. These vevs transform under the U(l) and 
R symmetries. As a consequence of the invariance, we find 

detM t 



es = f((h---,(m) + ^~ lo S 



+ ••• (9) 



M G detM^ 

where the first term represents the contributions from W sym breaking- It is only a function 
of U(l) and R invariant products of powers of vevs {Q}- The second term, coming from 
the Higgs sector, is discussed further in the next section and the ellipsis refers to the small 
additional contribution arising from the fermion mass sector. Note, M t , M' d only include 
those states, from 5s and 5s of SU(5) contained in W sym breaking and Wmggs, which mix with 
the Higgs sector (see the next section). 



B. The dependence of £3 on the Higgs sector 

We assume the Higgs sector includes any number of 10s with, for simplicity, only one 
of these, say 10i, coupling to light fermions. We also include the doublet-triplet splitting 
mechanism introduced by Dimopoulos-Wilczek |16|]. Accordingly, the terms in the super- 
space potential relevant for doublet - triplet splitting are of the form 

W d _ t = 5 2 K | (B - L)} 5i - §1 h ~ (B - L)\ 5 2 + ]T M v 5, 5, (10) 

1 1 i,j>2 

where a± | (B — L) is the vev of the field Ax in the 45 dimensional representation and 
(B — L) is the [Baryon - Lepton number] charge matrix (see eqn. (0)). Thus, since 
the doublets in 5i,5i C 10i have zero B — L, they remain massless. [Note, some of the 
5i , 5j states may come from 16 and 16 representations, respectively, in W sym breaking- We 
include, however, only those states which mix with 10i.] With this superpotential, we have 
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M' d = M[d], and the triplet mass matrix, M t , with non- vanishing determinant, includes the 
terms which mix states in 10i with 5j, 5j, for i, j > 2 and the sub-matrix M[t], where the 
matrix M[d], (M[£]) is given by M with Clebschs appropriate for doublets (triplets). 

Proton decay in this model is mediated by the Higgses in 10i. Hence, the coefficient of 
the resulting effective dimension 5 baryon violating operators ||[7| is given by the inverse of 
the triplet mass matrix in the 11 direction, i.e. 

_ detM [t] detM> 

where g = ^Mll • ^ e snow ; m appendix 2, that g = g(Ci, • ■ ■ , Cm) is a holomorphic function 
of the set of U(l) and R invariant products of powers of vevs. If we now define the effective 
triplet mass by M t 1 = (M t _1 )xi, we obtain the final form for the factor in eqn. (§) - 
detM* = * fl'(Ci) ' ' ' 5 Cm)- As a consequence, eqn. @ becomes 



6 3 = F(Cl,...,C m ) + ^l0g 



M t 



(11) 



M G 

where / and g are absorbed in the function F. We thus find that suppressing the pro- 
ton decay rate by increasing the ratio has the effect of increasing the value of e% and 
consequently increasing the value of a s .f\ 

C. An Example : £3 in the model of paper II 

As an example, for the particular model of paper II, there is only one independent 
invariant ratio of the GUT scale vevs given by ( = In principle, 63 could depend on an 
arbitrary function /(C)- However, when we put the effective determinants for this model in 
equation @ we find 

' Mt 
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3 This observation in the context of minimal SU(5), where in that case Mt is the color triplet Higgs 
mass, was discussed earlier by Hisano et al. jjj. 
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Remarkably, all dependence of 63 on the GUT scale vevs has dropped out except for the 
dependence on M t . Unfortunately, the large positive constant that appears in the expres- 
sion means that in order to get €3 negative, the proton's lifetime would be many orders 
of magnitude lower than the experimental bound. In addition, changing the Yukawa cou- 
pling coefficients of the terms of W sym breaking cannot cure this problem. For the most part, 
changing the Yukawa coupling constants of one of the terms of W sy mbreaking has the result 
of multiplying the effective determinants of the mass matrices M 7 in eqn. (^), for states 7 
in a complete SU(5) multiplet, all by the same amount. Thus, this multiplicative factor has 
no effect on 63. Hence, as a consequence of (0), we find that the symmetry breaking sector 
of the model of paper II is ruled out by the dual constraints coming from the low energy 
measurement of a s and the proton lifetime. 

Is it possible to find a symmetry breaking sector which has all the U(l) symmetries 
required for the naturalness of the theory and is consistent with the constraints of a s and r p ? 
In the next section we describe a new SO (10) SUSY GUT which satisfies all our criteria. In 
fact it may contain the minimal symmetry breaking sector consistent with the requirements 
of (1) obtaining the effective fermion mass operators of paper I, (2) "naturalness," and (3) 
retaining only those states at energies below Mq, which either have trivial SM charge or are 
contained in the minimal supersymmetric standard model [MSSM]. 

III. A COMPLETE SO(10) SUSY GUT 
A. The GUT symmetry breaking and Higgs sectors 

For this model, 

^^sym breaking 

±A[{Al + S 3 SA 1 + S 4 A 1 A 2 ) (13) 
+ A 2 (^ + S 1 A) + SA 2 
+ S'(SS 2 + A X A) + S 3 S' 2 
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where the fields {Ai, A 2 , A, A[}, {S, S'}, if), ip, {Si, ■ ■ ■ ,S^} are in the 45, 54, 16, 16 and 
1 representations, respectively. [Note, traces and contractions of indices are implicit.] The 
supersymmetric minimization condition = gives four discrete choices for the direction 
of the vev of A±. We will assume that nature chooses the B - L direction. The second term 
gives A a vev in the X direction, A 2 a vev in the Y direction, and ip and ip vevs in the 
right-handed neutrino-like direction (see eqn. fll4]) below). The third and fourth terms, and 
the S4 subterm of the first term, were added to give mass to all non-MSSM fields which are 
not in a singlet representation of the standard model gauge group. The last term was added 
in accordance with our "naturalness" criterion, namely that the theory should be the most 
general one consistent with the symmetries. The term S3S' 2 is consistent with the U(l) and 
R symmetries of the symmetry breaking sector of the theory that will be discussed below, 
and we are aware of no other additional symmetry that might exclude this term. Therefore 
the term must be included. 



The above vevs are given by 



(Ax) = ax 



r, = ai -(B-L) 



(14) 



(A 2 ) = a 2 
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V = a 2 -Y 
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(if,) = v |SU(5) singlet) 



(^) =tj |SU(5) singlet) 
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The vacuum minimization conditions are explicitly 

a\ = — sSs, S±a = -vv (15) 
2 

S1S2 + -a\a = 0, a2<Si + 2sa = 
5 

Using these equations, the set {a±, 02, a, v, v, S4} form a complete set of independent vari- 
ables. 

Two caveats - 

• Note at tree level the GUT symmetry breaking vevs are undetermined since the po- 
tential in these directions are both supersymmetric and flat. We will not discuss the process 
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for fixing these vevs in this paper. That analysis must necessarily include supersymmetry 
breaking effects as well as supergravity and radiative corrections. 

• We do not describe the source of supersymmetry breaking in this paper. Soft SUSY 
breaking operators are included at the GUT scale and renormalized down to the weak scale 
when making any comparison with the low energy data. 

The same doublet-triplet splitting mechanism that was used in paper II can be used in 
this model. Accordingly, the Higgs sector of the Lagrangian is given by 

Lm 99 s = [1(W0 2 + S 5 10j\ F + ^*10?b (16) 

where the first term is Wniggs and the second generates a \x term when the F component of 
the hidden sector field z gets a vev of order \Q 10 GeV. We then obtain \x = (F z )/M. 



1. U(l) and R symmetries of W symhreaking and e 3 

Wgymbreaking has a [U(l)] 4 x R symmetry as is summarized in table [TTI|. Up to arbi- 
trary Yukawa coupling coefficients assumed to be of 0(1), W sym breaking is the most general 
superspace potential consistent with these symmetries. Using ai,a,2,a,v,v, and <S 4 as inde- 
pendent variables, the only invariant under a [U(l)] xR rotation of the vevs is £ = -^i- 
After evaluating e% explicitly using equation (H) we obtain 
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(i-C) 

Before we can check whether the experimental measurement of a s (Mz) in this model is 
consistent with the non-observation of proton decay, we must discuss the fermion mass sector 
of the theory. The proton lifetime and branching ratios depend crucially on the couplings 
of the color triplet Higgses to fermions. In a theory of fermion masses, these couplings are 
related to the Yukawa couplings of fermions to Higgs doublets, but they are not identical. 
For example, the dimension 5 operators responsible for proton decay are given by 

i Q|Cq Q Q QC ql L (18) 
jjjT uC ue e uC ud d 
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where Cqq, Cql, C ue and C u d are 3x3 complex matrices. The matrices Cqq and g u , the up 
quark Yukawa matrix, contain the same independent parameters but the SO(10) Clebschs 
are different. An example is presented in the next section. 



B. Fermion mass sector 



We take the flavor sector of our model to be that of model 4 in paper I. Preliminary 



results of Blazek, Carena, Wagner and one of us (S.R.) [0 suggest that this model provides 
the best fit to the low energy data. 

The flavor sector is specified by a particular set of four operators {033,032, 022,012}. 
Three of these operators - O33, O32, and O12 - are uniquely specified by choosing model 
4. On the other hand there are 6 choices for operator O22, labelled (a, • • • , /), as all give 
identical entries in the charged fermion Yukawa matrices. As a result we can construct 6 
possible models 4(a, •••,/). 

The four effective fermion mass operators for model 4c are given by 

33 = 16 3 10! 16 3 (19) 

23 = 16 2 f 10i ^16 3 
22 = 16 2 £ 10, £ 16 2 

The superspace potential for the complete theory above the GUT scale which reproduces 
model 4c is given by (see fig. 1) 

(20) 

I63IO1I63 + \f1Al63 + ^ 1 A9 1 + il0i# 2 

+ f 2 i* 2 + \IM2I62 + \IM1I62 

9 

+ ^ 3 10l^4+ S M ^2(^^a) 

a=3 
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+ # 6 10i^ 7 + ^7^4^8 + ^8-4^9 + ^9^4l6i 



This superpotential is consistent with the symmetries discussed previously with the ad- 
dition of one new U(l) given in table |TTT| . However it is not the most general fermion sector 
consistent with these symmetries. In fact one and only one new operator must be added 

f 6 A 2 16 3 . (21) 

It is easy to see that this operator leads to one new effective operator at the GUT scale 
when heavy states are integrated out (see fig. 2). The new operator is 

°i3 = i6i(^) 3 10^) 163. (22) 

The complete model 4c is thus defined with this new operator and includes the operators in 
eqn. ([19]) and (E^). 



A fit of model 4c to the low energy data |T2| for certain ranges of soft SUSY breaking 
parameters agrees to better than la for all observables. A global \ 2 analysis of models 4, 
6 (a - f) (paper I) for all regions of parameter space is now underway fli~2"f . It is already 
clear, however, that an additional operator such as the 13 operator in eqn. ( p2[ ) is absolutely 
necessary to fit the data. Whether this model fits the data better than any other remains to 
be seen. For example, a different choice of 22 operator results in a different U(l) symmetry 
and thus by "naturalness," a distinct theory. In particular, we have checked that for model 
4b there are no new effective fermion mass operators generated, while for model 4a the new 
13 operator 

°- = 16 1 (£) 3 10i(^)l6 3 (23) 

is needed. 

Finally, consider the matrices relevant for proton decay for the case of model 4c. In 



particular, the matrix Cqq is given by 
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(24) 



This should be compared with the up quark Yukawa matrix in the same model given by 



9v 



' C \De iS \ 



C -\B 



(25) 



where the 7 Yukawa parameters A, B, C, D, E, <ft and 5 are obtained by fitting to quark 



and lepton masses and the CKM mixing angles [fL2]| . Note that the Clebschs in the two 
matrices differ by as much as a factor of 4. These Clebschs affect proton decay branching 
ratios. It is thus important to calculate the branching ratios in models which are consistent 
with the observed fermion masses and mixing angles. 

1. Additional threshold corrections at Mq 

The dominant effective operators at the scale Mq are obtained by integrating out states 
with mass greater than Mq- Higher order corrections to these operators are also obtained. 
These typically lead to 0(10%) corrections to the leading terms in the Yukawa matrices. Of 
course the terms in the Yukawa matrices will also receive corrections at one loop. We have 
neglected these corrections in the following analysis. 

2. Neutrino masses 
For completeness we include a minimal neutrino mass sector 

3 

W neutrino = ^16, Ni (26) 

i=l 

where the states N, i = l,---,3 are SO(10) singlets. This term has the effect of giving 
GUT scale Dirac masses to the right-handed neutrinos in the 16 's. Thus, the SM left-handed 
neutrinos are absolutely massless in this theory. 

If necessary, left-handed neutrinos can be given masses as described in paper II fllP 



However, in order to do so we must either break one linear combination of the U(l)s or 
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introduce additional SO (10) singlets. In either case, we must then check for "naturalness" 
and add any operator allowed by the symmetries of the theory. 



C. Symmetries 

This theory has 5 global U(l) symmetries and a global continuous R symmetry. The 



charges of most of the states are given in table The charges of the other states can easily 



be derived. We have checked our theory for "naturalness." We find that only 3 additional 
operators need to be added to the superpotential - 

<MiV>i<S2, tM'i16 3 S 3 , foAifoSs. 

These three operators have no direct effect on any observable properties since the vev of 
A[ vanishes. With the inclusion of these three operators the total superspace potential for 
model 4c is "natural" (i.e. no additional operators consistent with these symmetries are 
allowed) for all possible powers of the fields. In addition the theory has a matter reflection 
symmetry (see Dimopoulos and Georgi, ref. [1]]) which forbids dimension 4 baryon or lepton 
number violating operators. 

Some problems, however, remain to be solved in our model. Given the states and sym- 
metries, we find that f a @, the coefficient of the general gauge kinetic term, is trivial in this 
model. Thus we have not explicitly included the sector of the theory which generates gaug- 
ino masses once SUSY is broken. In addition, the U(l) symmetries of the theory are not 
sufficient to significantly constrain the Kahler potential. For example, terms such as 

jf 2 K SI, A 16a 

are allowed which mix light generations with heavy states. This term (and others like it) is 
allowed since it already appears in the structure of the Feynman diagrams of fig. 1. These off 
diagonal terms in the Kahler potential can affect fermion mass operators as well as introduce 
flavor changing neutral current processes at low energies [ T7fl . When deriving the effective 
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theory at Mq we have implicitly assumed that the Kahler potential is universal for all 16s 
in the theory. 

IV. RESULTS FOR a s (M z ) VS. PROTON DECAY IN MODEL 4C 

We now check whether the new model is consistent with a s and the proton lifetime. We 
have calculated the proton decay rate for the three dominant modes — K + u, tt + u, K°h + 
where for neutrino modes we sum over the three neutrino species. We have included both 
gluino and chargino loops, as well as LLLL, LLRR and RRRR operators in our analysis, 
where L(R) refers to left-handed (right-handed) fermion fields. We used values for the 
dimensionless Yukawa parameters at Mq which give predictions for fermion masses and 
mixing angles in excellent agreement with the data [0 . The values for soft SUSY breaking 
parameters are also consistent with electroweak symmetry breaking and the experimental 
measurement of b —>■ s + 7. We renormalized the dimensionless (dimensionful) parameters 
at two (one) loops to low energies in order to make contact with the data. A detailed 
report on this work is in preparation ||13|| . For the calculation presented here we take the 
effective Higgs triplet mass, M t = a\/S h = 4 x 10 19 GeV with a x = M G « 2 x 10 w GeV and 
1S5 = 10 13 GeV. This corresponds to light Higgs doublets with mass 10 13 GeV. 

Is it natural to have such light Higgs doublets? Are we populating the GUT desert? 
To address this question we should compare the Higgs doublet mass with the spectrum of 
masses for states in the symmetry breaking sector of the theory. These in fact range from 
10 13 - 10 20 GeV. Thus we have taken the doublet mass to lie at the lower bound of this 
GUT scale spectrum. This seems to be the only natural criteria for setting an upper bound 
on M t . 

The results for proton decay are given in table [IV] for two different values of soft SUSY 
breaking parameters and dimensionless couplings ||12|| . For comparison we also present the 
branching ratios obtained in a generic minimal SU(5) SUSY GUT from Hisano et al. 0. The 
soft SUSY breaking parameters at the GUT scale are for (case A) Mi/ 2 = 250, m = 750, ji = 
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64,m Hu = 904,m Hd = 1200 GeV and (case B) M 1/2 = 100, m = 3000, fi = 322, m Hu = 
4200, vn,n d = 3150 GeV with for both cases, A = and tan f3 = 53. The experimental lower 
bound on the proton lifetime into the mode K + v is 10 32 years |TJ|. Thus these values are 
consistent with the non-observation of proton decay to date. However, these values are to be 
considered as upper limits on the proton lifetime. In particular, r p Br~ 1 {p — > K + i?) scales 
as ( lo^GeV ' 003 /^ 3 ) 2 ' where Z 3 ' with values in the accepted range (3 = (0.003 - 0.03)GeU 3 



|19| , is a measure of the matrix element of a 3 quark operator between the proton state and 



the vacuum. 

Are these upper bounds on proton decay consistent with the experimental measurement 
of a a {Mz)l We have evaluated the expression for e 3 including the additional states with 
GUT scale masses contained in the fermion mass sector of the theory. These typically shift 
the value of tz by a small amount in the positive direction. We find that for typical values 
of ax, a2, and S4 around the GUT scale, we can obtain 63 negative for M t of order 10 19 GeV. 
For example, with d\ = 2a 2 = 2S4 = | = Mq we have £ = 16 and €3 ~ —0.030, which 
includes a positive contribution of 0.005 from the fermion mass sector. 

We have also checked that for these values of the parameters the gauge coupling satisfies 
one loop perturbativity up to M = 10 18 GeV with a G (M) = 0.39. Note that above the 
GUT scale we use the threshold boundary condition 

olg\M g ) = a5 1 +A 2 (M G ) 



V. CONCLUSION 

We have presented a new complete SO(10) SUSY GUT. This theory has several inter- 
esting features. The superspace potential is "natural" to all orders in the fields. It contains 
what may possibly be the minimal GUT symmetry breaking sector necessary to obtain the 
desired adjoint vevs consistent with (1) "naturalness" and (2) fermion masses. 
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We have also shown that one loop GUT scale threshold corrections are a significant 
constraint on the GUT symmetry breaking sector of the theory. This constraint, for example, 
is sufficient to rule out the model of paper II ||10| . The one loop threshold corrections relate 
the prediction for a s (Mz) to the proton lifetime. 

Finally we have calculated one loop GUT scale threshold corrections to gauge couplings 
and proton decay rates to different final states in the new model. The results are consistent 
with the low energy measurement of a s (Mz)- The proton decay branching ratios provide a 
powerful test for theories of fermion masses. There is reasonable hope that new results from 
SuperKamiokande or Icarus could confirm or rule out these theories. 

It would be presumptious of us to conclude without discussing some of the open ques- 
tions. We have not discussed the origin of SUSY breaking nor how it feeds into the visible 
sector with the one exception of the \x term which we have nominally considered. We also 
do not discuss what determines the GUT vevs. At tree level, in the globally supersymmet- 
ric theory considered here and neglecting SUSY breaking, these are flat directions of the 
potential. Finally, and perhaps most seriously, the symmetries discussed in this paper do 
not significantly constrain the Kahler potential. Flavor mixing in the Kahler potential could 
lead to dangerous flavor changing neutral current processes. We have assumed the trivial 
universal Kahler potential in our analysis. 
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APPENDIX 1 



Proof of equation (|6]) 

We define the GUT scale M G as the point where a G = oi\(M G ) = a^Ma)- This means 
that Ai(M G ) = A 2 (M G ). We then define the relative shift in a 3 (M G ) by 

e 3 = (a 3 (M G ) - a G )/a G (27) 



a 3 (M G ) (A 3 - A 



l\M G , 



We then have 

A 3 - Ail^ = ^-E(62 " &?) log |det'M 7 | - ]T(^ - 6<»gM G ) (28) 

27F 7 7 

where n 7 = the mass dimension of det'M 7 and det'M 7 is defined in eqn. (|7|), except for 
det'Ma where it will be convenient to define Md by det Md = M G detM^, and det'M^ = 
detMrf. This redefinition does not affect eqn. (|6|). Note, the matrix Md is defined such that 
n t = n d . 

In addition, 

Ai|m g = A 2 |m g 

which implies 

J2(bJ - bl) log |det'M 7 | = (£(bj - bl)n y ) log M G (29) 

7 7 

Substituting for log M G we obtain 

^3 - fpffi^ - bl) log |det'M 7 | - |^|^ £(*Z - bl) log |det'M 7 |} 

= ^— E(^ c 23 + b\c zl + bjc 13 ) log |det'M 7 | (30) 

where cy = I] 7 (&7 — &J)n 7 . To evaluate the c^s, define ns4, n 45 and n w to be the number of 
54, 45, and 10 representations in the theory, respectively, and uiq and njg to be the number 
of supermassive 16 and 16 representations, respectively. For any SO(10) model built with 
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only 1, 10, 16, 16, 45, and 54 representations and one pair of light Higgs doublets, we have 









n g 




= n 45 + n 54 




n x 


= n 45 + n 54 - 3 


n u 


= n e 


= n Ab + nie - 3 


n s 


= n a 


= ™54 




n q 


= n 45 + n 54 + n 16 


n d 


= n t 


= nio + n 46 



Evaluating C23 explicitly 

c 2 3 = (4n s + 3n q + 2n w + 3n x + n d ) - (5n s + 2n q + n u + 3ri 9 + 2n a . + n t ) 
= {4(n 54 ) + 3(n 54 + n 45 + n w - 3) + 2(n 54 + ra 45 ) + 
3(n 54 + n 45 - 3) + (ni + n 16 ) 

(31) 

-[5(n 54 ) + 2(n 54 + n 45 + rii 6 - 3) + (n 45 + n i6 - 3) + 3(n 54 + n 45 ) + 
2(n 54 + n 45 - 3) + (n w + n w )]} 
= -3 

Similarly 

c 3 i = 9 (32) 
C12 = -6 

Thus, the CijS are completely independent of the number of fields in any theory built with 
Is, 10s, 16s, 16s, 45s, and 54s. Plugging eqns. ( |5T] ) and ( |3"2"D into eqn. Q3"0|), eqn. (|^) readily 
follows. 

APPENDIX 2 

U(l) symmetries and the dependence of e 3 on GUT scale vevs 

In this appendix we prove that the contribution to £3 from the GUT symmetry breaking 
sector is only a function of U(l) and R invariant products of powers of vevs. The proof relies 
on two facts: 
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1. the effective determinants of mass matrices are holomorphic functions of the symmetry 
breaking vevs, and 

2. the effective determinants have simple phase rotations under U(l) and R symmetry 
transformations. 

Note, since the effective determinants are independent of the conjugates of vevs, the U(l) 
and R invariance of e 3 is very restrictive. 

We first discuss the case for mass matrices which do not include vector multiplets, fol- 
lowed immediately by the case for mass matrices including vector multiplets. Note, in the 
first case the mass matrices themselves are, by construction, holomorphic functions of vevs. 
This is however not true for the latter case which is why it requires a separate discussion. 

Consider a general superspace potential W($i, <E> 2 , • • • , &n) whose superfields rotate un- 
der a U(l) symmetry as 

By defining the shifted fields = — and expanding the superspace potential about 
($) we can find the fermion mass matrices. 

Wfa + (<£>!>, ...) = ... + £ VW($i>, (*2>, • • 0^7 + • • • (33) 

7 

Now consider what would happen if the superfields received a different vacuum expectation 
value, 

($j) new = e iQ 3 e ($ j ) old 
Under this change, the mass matrices would change. 

Wfa + e iQl6 '<$!), ...) = ... + J2^rn% + ... (34) 

7 

where 

m* = m 7 (e < « lfl ($ 1 ),e i ^($ 2 >,...) 

However, if we rotate the shifted fields $j by e 1 ® 1 ®, the superspace potential will be invariant 
under the combined rotation of $ and (<&). 
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W(e iQld ^ l + e iQl9 ($i), . . . , e iQN6 $ N + e iQNd (<f> N )) 



= . . . + E 7 ^ 7 



V 



V 



v> 7 + . . . 



(35) 



= . . . + J2-y ^7 m 7 ^7 + • • • 

Therefore, m® rotates in a very simple way. 

/ e -iQi6 \ ( e ~ iQl 



m" 



to- 



Thus, 



detm^J = e" 2if? ^ Q detm 7 



(36) 



(37) 



where the sum is over all fields that have columns in the mass matrix. 

Similar arguments can show that under an R symmetry rotation, ($) — > e tQd (§), to 7 is 
equal to 



I e ~iQit 



V 



\ 



to- 



V 



-iyic 



(38) 



where is the charge of the superspace potential under the R symmetry. Therefore, 



detrn^ = jQwrn-ne^Q de tm 7 



(39) 



where N = dimm. 



The situation is a bit more complicated for mass matrices which receive contributions 
from vector multiplets. The proof that the determinants and hence the effective determinants 
have simple phase rotations under the U(l) and R symmetry transformations can readily be 
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extended to these mass matrices. However, since the entries in the gaugino-chiral fermion 
mixing rows are actually the complex conjugates of vevs, the determinants of these matrices 
are not holomorphic. However, in the following we prove that the effective determinants of 
these mass matrices, which include vector multiplets, are in fact holomorphic functions of 
the vevs. 

Since the would-be Goldstone fermion states are perpendicular to the massive chiral 
fermion states, any mass matrix containing gaugino-chiral fermion mixing can be written in 
the form 

/ x\ x* 2 xl 

X\ Cn C12 C13 
X* 2 C 2 i C12 c 23 
^3 c 31 c 12 c 33 



X* N \ 
ClN 

C3N 



(40) 



\ X*N C N1 CN2 C7V3 " " " C NN J 

with J2jCijXj = for all i, J2iCijXi = for all j, and the c^s, xs, and xs are functions of 
the vevs but not of their conjugates. The rows and columns of the mass matrix can be 
rearranged so that xn and x N are not zero. However, the determinant of this matrix can be 
reduced by elementary row and column operations. Namely, by adding to the last column 
the second column multiplied by ^ plus the third column multiplied by and so forth, 
the determinant of the mass matrix becomes 

n -r* <r* • • • -r* t* -\- ^1^1 + 4- ^- lXiv ' 1 

u x l x 2 x 7V _ 1 x N -t- xn -t- . . . -+- xn 



x\ cn C12 • • • Ci )JV _i 

^2 c 21 C12 • • • C 2iA r_i 

Xl C31 C12 • • • C 3iA r_i 

^AT Cati Cat 2 • • • CAf 5 AT_i 
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Doing the analogous operation on the rows, the determinant becomes 






', 




X JV-1 


1 v^W 




Cll 


C12 


Cl,AT-l 





x 2 


C21 


C12 


C2,JV-1 






X JV-1 Cjv-1,1 Cjv-1,2 

=^Ef^x fc 



Thus, the determinant is equal to 



x k X k Y,k x *k x k 



X N 



X N 



CJV— 1,1 



CN-1,N-1 



Cl,JV-l 
■ C/V-1JV-1 








(41) 



N N 

= (J2 x k x k)(Y,x*k x k)f(vevs) 
k k 

where / is a holomorphic function of the vevs. By setting v = v, Ti will equal Xj for all i 



and the mass of the vector multiplet is yJ2k x \ x k- Therefore, the determinant is equal 
to M^ ector mult iplet times / and the effective determinant is just the function /. Note, 
M vector muitipiet is thus always canceled from the denominator of the effective determinant 
(eqn. ||) and no conjugates of vevs can appear in the effective determinant. Thus, the 
effective determinants for mass matrices containing gauginos are holomorphic functions of 
the vevs; just as those discussed earlier for mass matrices which do not have gaugino-chiral 
fermion mixing entries. 

These simple transformation properties of the mass matrices, under U(l) rotations of 
the vevs, have significant consequences for the form of £3. Consider the following expression 
entering 63 (see eqn. H) — 



-logdet'— 2- 

2 & Mq 



3. ,M tt 
- log det 

2 & Mq 



33, 

H log det 

10 5 



,M a 21 M a 

log det 

M G 10 M G 



(42) 



+ 10 1Ogdet Mc + 10 1Ogdet Mc 



6, , ,M q 
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It is now easy to show that it is invariant under the U(l) and R symmetries. Namely, the 
U(l) charge of the determinant of M w (eqn. |3?|) is equal to the charge of the determinant of 
M g which is equal to -2 times the sum of U(l) charges of all fields in the 24 representation 
of SU(5). Therefore, the U(l) rotation of detM s will cancel the rotation of detM w in ex- 
pression, eqn. (£E^). In addition, we note that the U(l) charges of the effective determinants 
of {M s , M a }, {M u , M e }, and M q equal -1 times the sum of U(l) charges of all fields in 
the {15 and 15}; {10 and 10}, and {10, 15, 10 and 15} representations of SU(5), respec- 
tively. Therefore, the U(l) rotation of det'M g is canceled by the rotations of the effective 
determinants of M u , M e , M s , and M a in expression, eqn. (f42|) . Similar arguments show that 



the expression in eqn. (|^) is invariant under an R symmetry rotation. Thus finally we 



arrive at the conclusion that the expression in eqn. (42) is invariant under the U(l) and R 
symmetries of W sym .breaking for any superspace potential built with 1, 10, 16, 16, 45, and 
54 representations of SO (10). Moreover, since expression, eqn. (|42"D , is holomorphic, the 
contribution of the GUT symmetry breaking sector to 63 is only a function of U(l) and R 
invariant products of powers of vevs. 

The same is not true for the contributions from either the Higgs or fermion mass sectors. 
This is because both the Higgs and fermion mass sectors contain massless states that must 
be projected out of the mass matrices before the effective determinants are taken. After 
this projection, the determinants of the resulting mass matrices are no longer holomorphic 
functions of the vevs. Nevertheless for the Higgs sector we can prove a similar but limited 
result, namely, that g = ^^uj = <7(Ci> ' ' ' > Cm); i- e -> 9 is a function of U(l) invariants only. 
By eqns. ( j37[) and fl39|), we see that detM[i] and detM[d] transform in the same way under 
the U(l) and R symmetries; therefore the ratio is invariant. 
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TABLES 



TABLE I. SU(5) decomposition of SO(10) fields in our theory 

54 -> 24 + 15 + T5 
45 -> 24 + 10 + TO + 1 
16 -> 10 + 5 + 1 
10^5 + 5 



TABLE H. Notation used for states in different charge sectors 



SU(3) x SU(2) x U(l) representation 


name 


appears in SU(5) representation 


(8,1,0) 


9 


24 


(1,3,0) 


w 


24 


(3, 2, -|); (3,2,|) 


X] X 


24 


(3,i,|); (3,i,-|) 


u; u 


TO; 10 


(1,1,-2); (1,1,2) 


e; e 


TO; 10 


(3,2, |); (3,2,-|) 


q;q 


10,15; T0,T5 


(6,1,-|); (6,1,|) 


s; s 


15;T5 


(1,3,2); (1,3,-2) 


a, a 


15;T5 


(3,1,-|); (3,1,|) 


t;t 


5;5 


(1,2,1); (1,2,-1) 


d] d 


5;5 
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TABLE III. U(l) and R charges of the new model 



field 


Till l rharffp 


field 


T J ( 1 I rharffe 




n n n n ni 

U, U, U, V) 


A' 


('_'? n 01 


An 
-ri-2 


fo 1 n n nl 

yv, i, u, u, 


A 


fn 1 n 01 

yv, u, ±, u, 




fn n -2 n rn 




(—] n — i 01 

l 1, U, ±, U, V) 


oh 


l U , ->-; U ; - 1 ) U J 


o~h 


(n n n — 1 nl 

l^u, U, U, 1,UJ 








n n ^ n r\\ 

[1, U, O, U, V) 


<->3 


l z ) U, Z, U, UJ 


9, 
O4 


(9 1 n n nl 

l z ) 1,U, U,UJ 


c 


i' 9 n n n a\ 

[Z, U, U, U, — 4J 


ID3 


fn n n n 11 

l U ; U, U, U, 1J 


16 2 


(-1,-1,2,0,1) 


16l 


(-2,-5,7,0,1) 


10i 


(0,0,0,0,-2) 


10 2 


(-1,0,0,0,2) 


z 


(0,0,0,0,-4) 


N 3 


(0,0,0,1,-1) 


N 2 


(1,1,-2,1,-1) 


Ni 


(2,5,-7,1,-1) 



All fields, except A[ and z, have R charge 1. A[ (z) has R charge (2). The superpotential has R charge 
3. The first 4 charges listed above contribute to states in the GUT symmetry breaking sector. 



TABLE IV. Proton decay results for model 4c 





A 


B 


C 


t p Br~ x {p -> iT + z/)/(10 3 Vs-) 
r(p^7i-+p) 

T{p^K+P) 
T{p-^K+u) 


66 
0.28 
0.0040 


3000 
1.18 
0.0045 


0.5 
0.0007 



For M t = 10 19 GeV,P = 0.003 GeV 3 . Note, t p Br'^p -» K+v) scales as ( ni^/' 003 ^ 3 ) 2 . The soft 
SUSY breaking parameters for case A and B are described in the text. In C, we show the results from 
Hisano et al. for comparison. Those results are from a generic minimal SU(5) SUSY GUT. 
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